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ABSTRACT 

We present coefficients for tlie calculation of the continuous emission spectra of H I, He I 
and He II due to electron-ion recombination. Coefficients are given for photon energies from 
the first ionization threshold for each ion to the n — 2Q threshold of hydrogen (36.5/i), and 
for temperatures 100 K <T< 10^ K. The emission coefficients for He I are derived from 
accurate ab initio photoionization data. The coefficients are scaled in such a way that they may 
be interpolated by a simple scheme with uncertainties less than 1 % in the whole temperature 
and wavelength domain. The data are suitable for incorporation into photoionisation/plasma 
codes and should aid with the interpretation of spectra from the very cold ionised gas phase 
inferred to exist in a number of gaseous clouds. 
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!_H 1 INTRODUCTION 



Analyses of optical recombination line (ORL) spectra of pho- 
toionised regions have suggested the existence of cold ionised gas 
(100 K-2000K), mixed within a warmer component at more typi- 
cal nebular electron temperatures (8000 K-10000 K). The pockets 
of cold ionised gas of high metal content are invoked as a possible 
explanation of the long-standing problem of the discrepancy be- 
tween elemental abundances derived from ORLs and those derived 
from collisionally excited lines (CELs). In this scenario, ORLs 
and CELs are preferentially emitted by the cold and warm phases, 
respectively (for a recent review see Liu, 2002); this problem is 
closely linked to the observation that Balmer jump temperatures of 
H II regions and Planetary Nebulae are systematically lower than 
those derived from the [O III] nebular to Auroral line ratio (Liu & 
Danziger 1993). 

Understanding the complex spectra arising from such regions 
relies on the construction of detailed photoionisation models able 
to account for all gas phases that may be present. Currently, one of 
the major limitations in such modelling is the lack of an accurate 
atomic data set extending to such low temperatures. The need of 
low-temperature effective recombination coefficient for the calcu- 
lation of recombination lines for metals is apparent. However, in 
addition to the discrete emission line spectrum, a continuous emis- 
sion is also produced by the ionised gas, mainly due to free-bound 
recombination processes of hydrogen and helium ions, free-free 
transitions in the Coulomb fields of H"'", He^ and He^"'' and two- 



photon decay of the 2 level of H I and He II, and, less im- 

portantly of the 2 ^So level of He I. Accurate continuous emission 
coefficients are essential for the correct prediction of the Balmer 
jump by photoionisation codes. 

The importance of the continuum processes listed above has 
long been known and emission coefficients have been tabulated (see 
e.g. Seaton 1955, 1960; Brown & Mathews 1970), for a range of 
temperatures and wavelengths mainly aimed at the study of op- 
tical data for classical H II regions. Ferland (1980) derived H I 
and He II continuous emission and recombination coefficients for a 
wider range of temperatures (500 K-2-10® K) and wavelengths to 
aid the interpretation of ultraviolet (UV) and infrared (IR) as well 
as optical data from nova ejecta. 

In this work we present new calculations of the H I, He I and 
He II continuous emission coefficients due to free-bound recom- 
bination, over temperatures ranging from 100 K to 10^ K, for the 
physical conditions thought likely to occur in chemically inhomo- 
geneous regions, which may include pockets of cold ionised mate- 
rial intermixed within typical nebular gas. 



2 CALCULATIONS METHOD 

Using the Saha-Boltzmann equation and the Milne relation, we 
may express the continuous emission coefficient ^{i^) correspond- 
ing to the recombination process 
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Figure 1. Frequency-dependent continuum emission coefficients 7(v) for H I {solid line), ^(y) for He I (dotted line) and ^(y) for He II (dashed line), 
calculated at T = 10000 K (left panel, to be compared with Figure 1 from Brown & Matthews, 1970) and T = 100 K (right panel). 



in terms of the photoionization cross-section av{X*) 

where w+ and ui* are the statistical weights of the recombining ion 
initial state and final state respectively, and e is the free electron 
energy. In terms of ^{v), the energy emitted per unit volume per 
unit time in frequency interval to + is N cN {X^)'^(i>)dv . 

For the contribution 7n(i^) at frequency ;/ from recombina- 
tions to states of principal quantum number n of H" and He"*" we 
compute the necessary energy-dependent photoionization cross- 
sections using the hydrogenic codes described by Storey & Hum- 
mer (1991). The total emission coefficient 7(1^) is then 

00 

no 

where in practice we truncate the sum at n = 200 for H I and 
n = 350 for He 11, which is sufficient to ensure convergence to all 
figures given in the Tables. 

For recombination to atomic helium we use the ab initio cal- 
culated photoionization cross-sections described by Hummer & 
Storey (1998). Cross-section data are available for n^S, n^S, n^P°, 
n^P°, n^D, n^D, n^F° and n^F° states with Z -f 1 < n < 20. For 
higher values of n and for Z > 3 we use hydrogenic data. Thus 
for n > 5 the emission coefficient has nine distinct thresholds for 
each n corresponding to the eight separate terms plus a threshold at 
the hydrogenic energy. Hummer & Storey (1998) showed that the 
results of their ab initio calculation of the photoionization cross- 
sections are in better agreement at threshold with the highly ac- 
curate bound-bound calculations of Drake (1996) than any of the 
other methods used to compute helium recombination processes. 

Contributions due to bremsstrahlung emission of a 
Maxwellian distribution of electrons in the Coulomb fields 
of hydrogen and helium ions are not included in our results, 
but can be obtained readily using (e.g.) equation 2 of Brown & 
Mathews (1970) or the free-free computer code published by 
Storey & Hummer (1991). 

Contributions to the continuum emission from two-photon 
emission are also not included but can be computed using the for- 
mulae of Nussbaumer & Schmutz (1984). 



3 RESULTS 

The full tables of continuous emission coefficients, Tables 3, 4 and 
5 for H I, He 11 and He I respectively are available in electronic 
form only. In Table 1 we show an extract from the table for H I. Co- 
efficients are tabulated for log T[K] = 2.0(0.1)5.0 and for photon 
energies (in Rydbergs) from just less than the ground state thresh- 
old energy to just above the energy of the threshold at n = 20. 
Values are tabulated on either side of each threshold and at some 
additional nodal points inserted to make interpolation sufficiently 
accurate. At low temperatures the coefficients fall rapidly and ex- 
ponentially at energies above each threshold making interpolation 
difficult so rather than ^(v) we tabulate 7^ {v) defined by 

7t(z.) = 7(,.) 10^4 T^/ V^/'^^ (1) 
= 7(,.) 10*° t^/^e^^-^««^^^«/* (2) 

where t = T[if]/10*, Ai5 is the difference between the photon 
energy, hi/, and the energy of the nearest threshold of lower en- 
ergy and /S.Er is the same energy in Rydberg units. Thresholds are 
indicated by index 1 and additional nodal points by index 0. 

The recommended procedure for deriving the emission coef- 
ficient at a given temperature and photon energy is to interpolate 
linearly in the appropriate table in the variables log T and pho- 
ton energy to obtain the scaled coefficient 7^(i^). Equation 1 can 
then be applied to obtain 7(1^). Table 2 gives exact values of 7(1^) 
and values derived from the recorrmiended interpolation scheme for 
each ion at a range of photon energies. 

Comparison of the tabulated values with calculations per- 
formed on a finer frequency grid and temperature grid show that 
linear interpolation in log T and photon energy yields accurate re- 
sults with maximvun deviations of 1% and average deviations much 
smaller than this. 

Figure 1 shows the non-scaled continuum emission coeffi- 
cients, 'yiy) for H I, He I and He II in the optical wavelength 
range for a temperature of 10000 K (left panel) and 100 K (right 
panel). The left panel of figure 1 is directly comparable to figure 1 
of Brown & Mathews (1970). We can only compare the magnitude 
of the discontinuity at each threshold directly with the results of 
Brown & Matthews (1970), since their results incorporate all con- 
tinuum processes while ours only deal with free-bound processes. 
Comparing results for H I and He 11 we find a maximum difference 



Table 1. [Abridged] Scaled H i free-bound continuum recombination coefficients, "/''{i') for H i [erg cm ^ sec ^ Hz 
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Table 2. Representative values of the exact and interpolated coefficients, 7(1/) 



V [Hz] 



E [Ryd] 



= [A] 



H i(exact) 



7(!/) [10-40 erg cm3 s"! Hz"!] 
Hl(interp) Hel(exact) Hel(interp) Hell(exact) Hell(interp) 



1.000(+14) 
2.000(+14) 
3.000(+14) 
4.000(+14) 
5.000(+14) 
6.000(+14) 
7.000(+14) 
8.000(+14) 
9.000(4-14) 
1.000(+15) 



0.030397 
0.060793 
0.091190 
0.121586 
0.151983 
0.182380 
0.212776 
0.243173 
0.273569 
0.303966 



29979.246 
14989.623 
9993.082 
7494.811 
5995.849 
4996.541 
4282.749 
3747.406 
3331.027 
2997.925 



2.415(+01) 
2.852(-|-00) 
1.653(+00) 
4.442(+01) 
6.996(-01) 
1.093(-02) 
1.698(-04) 
2.629(-06) 
2.376(+01) 
3.704(-01) 



2.418(+01) 
2.858(-|-00) 
1.655(+00) 
4.442(+01) 
6.988(-01) 
1.092(-02) 
1.694(-04) 
2.628(-06) 
2.376(+00) 
3.702(-01) 



2.408(+01) 
3.026(-|-00) 
2.067(+00) 
5.948(+01) 
2.848(-l-00) 
5.478(-02) 
1.015(-03) 
1.828(-05) 
1.795(+02) 
1.097(+01) 



2.408(+01) 
3.028(-|-00) 
2.068(+00) 
5.947(+01) 
2.856(-l-00) 
5.495(-02) 
1.017(-03) 
1.831(-05) 
1.795(+02) 
1.100(+01 



7.393(+01) 
3.222(-t-01) 
7.034(+01) 
9.356(+01) 
1.457(+00) 
3.098(+01) 
4.807(-01) 
7.437(-03) 
5.054(+01) 
7.823(-01) 



7.412(+01) 
3.233(-t-01) 
7.047(+01) 
9.366(+01) 
1.458(+00) 
3.098(-|-01) 
4.805(-01) 
7.438(-03) 
5.054(+01) 
7.821(-01) 



of 1% for H I at the Balmer threshold and at the lowest tempera- 
ture tabulated by them of 4000K. The difference is attributable to 
Brown & Matthews using an approximate expression for the hy- 
drogenic threshold photoionization cross-sections while we use the 
exact expressions incorporated in the codes of Storey & Hummer 
(1991). A similar comparison for He I shows that for all thresholds 
except that corresponding to the 3 state at an air wavelength 
of 7849A, our results differ by no more than 2.1% from those of 
Brown & Matthews (1970). The differences that do exist are due to 
the approximate method used by Brown & Matthews to calculate 
the helium photoionization cross-sections. In the case of the 3 "^P" 
threshold, we find much larger differences reaching 21% at 4000K. 
This is almost certainly due to a numerical error in the work of 
Brown & Matthews, since the magnitude of the discontinuity at the 
3 threshold does not obey the correct scaling with temperature 
in their work. 

We also compared our values of 'yiy) for H I and He II with 
those pubUshed by Ferland (1980) and found, in general, good 
agreement, with typical deviations of the order of 2-5% in the over- 
lapping temperature range. 

The electronic tables are structured as follows: node/threshold 
(0/1) indices are given in column 1, the photon energies in Ryd 
are given in column 2, the scaled free-bound emission coefficients 
for temperatures in the lOOK-lOOOOOK range are given in the re- 
maining columns. For He I and He II coefficients are tabulated for 
log T[K] = 2.0(0.1)5.0 while for He I, log r[K] = 2.0(0.04)5.0. 
The maximum photon energy is slightly less than the ground state 
ionization energy for each ion and the minimum photon energy cor- 
responds to the n = 20 threshold in H I. 



data are presented under an interpolation scheme that allows esti- 
mates to be obtained over the entire temperature and energy range 
presented with less than 1% uncertainty. 
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4 CONCLUSIONS 

We have presented new calculations of the free-bound continuous 
emission coefficients for the hydrogen and helium ions. The results 
for He I are derived from accurate ab initio photoionization cross- 
section data. The coefficients are given for a wide range of tempera- 
tures and frequencies, extending to the previously unexplored very 
low-temperatures regime. This is needed for the interpretation of 
spectra from very cold ionised gas, which has been inferred to ex- 
ist, possibly in the form of density/chemical inhomogeneities, from 
observations of ORLs in H 11 regions and planetary nebulae. The 



